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ABSTRACT 


stuJv  is  made  of  a  set  of  wel  1- recorded  Aleutian  Islands 
earthquakes  and  the  nuclear  explosion  LONG  SHOT.  Whereas  past 
studies'^'  travel- t ine  anomalies  and  location  errors  have  been 
based  on  seismic  signals  from  explosions,  the  main  object  of 
the  current  study  has  been  to  extend  our  working  principles  to 
seismic  data  from  earthquakes.  Using  a  selected  teleseisnic 
station  network,  travel-time  anomalies  were  computed  and  various 
techniques  applied  to  achieve  location  consistency  across  the 
entire  region.  Basically,  the  techniques  involve  deriving  space 
functions  for  the  observed  station  anomalies  and  imposing  the 
following  criteria  for  success:  tight  clusters  of  locations  for 
each  event  made  with  network  subsets,  acceptably  low  standard 
deviations  from  the  least-squares  solutions,  and  reasonable 
station  anomaly  functions. 

Using  LONG  SHOT  as  a  bias  control  point,  a  set  of  segmented 
constant  anomalies  were  derived  for  the  network  of  stations  used. 
Neglecting  depth,  the  relative  location  accuracy  is  believed  to 
be  about  5  km;  standard  deviations  of  solutions  generally  were 
reduced  to  less  than  0.25  sec,  a  value  accepted  as  being  due  to 
reading  and  timing  error.  Compared  with  published  locations,  the 
events  in  the  Rat  Islands  and  Near  Islands  as  ?  result  of  applying 
the  technique  shift  20-50  km  southerly,  and  those  in  the  Fox  Islands 
and  Andreanof  Islands  shift  10-20  km  southerly. 

To  demonstrate  the  difficulties  of  locating  events  in  the 
Aleutian  Islands,  well-distributed  sub-networks  of  a  529-station 
network  were  used  to  locate  LONG  SHOT.  The  locations  obtained, 
any  of  which  would  have  been  accepted  as  not  unusual  because  of 
network  coverage,  exhibit  errors  as  large  as  45  km  with  two- 
quadrant  networks  and  160  km  with  single-quadrant  networks;  travel 


tine  anonalies  computed  from  the  various  solutions  have  a  range 
of  more  than  11  sec  at  station  N'P-NT. 

To  verify  the  techniques  used,  a  study  is  also  made  of  a 
set  of  hypothetical  events  and  station  anomaly  functions.  It  is 
shown  that  (lj  even  though  a  constant  network  is  used  and  a 
constant  anomaly  is  the  only  error,  relative  accuracy  may  suffer 
due  to  a  network  effect  caused  by  the  nonlinear  travel  -  tine/di  stance 
relationships;  and  (2)  the  techniques  of  funct ional i c i ng  the  station 
anonalies  are  valid  and  could"be  implemented  as  recommended  for 
studies  in  selected  regions. 
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INTRODUCTION 


The  main  purpose  of  the  present  study  is  to  apply  to  earth¬ 
quakes  those  techniques  known  to  be  valid  for  accurately  locating 
underground,  explosions  within  a  snail  region  using  travel-tine 
anomalies  and  teleseisnic  networks.  Additionally,  we  hope  to 
demonstrate  techniques  for  extending  the  region  over  which 
functional  station  anomalies  can  be  derived  so  as  to  provide 
more  accurate  locations  than  those  obtained  fron  simple  least- 
squaies  schemes. 

Heretofore,  most  seisnic  studies  dealing  with  travel -tines , 
location  accuracies,  and  "station  correction^"  have  used  explosions 
to  provide  the  event  data  set,  because  only  then  can  one  be  certain 
of  the  hypocenter  and,  hence,  of  actual  travel  tines,  origin 
tines,  station  anomalies,  and  the  effects  of  an  inadequate  earth 
model.  Such  studies  could  meaningfully  investigate  tine  deviations 


on  the  order  of  0.2-0. 3  sec.  ShMl&vone  include  earthquakes  in  the 
data  set,  conclusions  regarding  location  accuracy  and  travel-tine 
anomaly  stability  froiii  one  region  to  another  always  remain  in 
question.  This  stems  fron  the  fact  that  errors  in  the  location  of 
an  event  assumed  correct  can  produce  errors  10-20  tines  the  site 
of  the  tine  anomaly  or  of  the  anomaly  change  one, is  trying  to 
observe  (Chiburis  and  Dean,  1965;  Chiburis ,  1966). 


Earthquakes  have  been  used  generally  to  obtain  first-order 
.effects  of  inhomogeneity,  principally  in  the  vertical  velocity 
distribution  of  the  earth  compared  to  a  standard  earth  model.  Two 
such  models  are  presently  in  common  use:  the  Herrin  (1968)  model 
akd  the  classical  Jef f revs-Bullen  model.  A  nunber  of  other  models 
have  been  proposed  (e.g.,  Archambeau  et  al.,  1969)  to  provide 

finer  details  in  the  velocity  distribution  usually  enphasiring 

/ 


( 


one  or  core  low-velocity  layers  or  transition  zones  not  included 
in  the  Herrin  and  Jeffreys-Bullen  models. 

It  is  now  clear  that  even  these  improved  models  are  inadequate 
to  give  tine-distance  relations  sufficiently  accurate  for  locating 
large  events  nuch  better  than  20-25  km  on  the  average.  For  small 
events,  and  correspondingly  larger  reading  errors,  the  location 
errors  can  easily  be  50  kn  or  nore. 

Among  others,  Chiburis  and  Ahner  (1970)  pointed  out,  that 
the  probable  causes  of  the  time  anomalies  are  due  to  unknown 
lateral  inhonogenei ties  in  the  mantle  (all  the  way  down  to  the 
core)  and  to  complexities  in  the  upper  mantle  and  crust  in  the 
vicinity  of  the  recording  stations.  It  is  doubtful  that  we  will 
ever  be  able  to  devise  a  three-dimensional  world  model  accurate 
enough  to  account  for  observed  anomalies.  However,  in  special 
cases  one  does  not  need  a  completely  accurate  model  simply  in 
order  to  locate  events:  anomalies  can  be  determined  from  an 
explosion  with  known  location  and  applied  to  subsequently  recorded 
explosions  in  the  same  region.  In  this  way,  locations  can  be 
obtained  with  an  absolute  accuracy  of  2-3  km  (Chiburis,  1968).  In 
fact,  the  exact  location  of  the  calibration  event  is  .immaterial; 
the  accuracy  o£  other  events  relative  to  it  remains  the  same 
(Chiburis  and  .^hner,  1970).  As  a  result  of  doing  this  for  many 
regions,  valuable  information  about  the  spatial  patterns  of  event 
occurrence  and  the  distribution  of  energy  release  could  be 
provided.  Most  important,  perhaps,  is  that  location  consistency 
among  many  events  could  be  achieved,  region  by  region,  with  the 
result  that  fewer  events  would  have  to  be  specially  analyzed- or 
would  have  to  be  excluded  from  bulletins  simply  because  of  poor 
time  fits;  furthermore,  those  events  included  would  display 
significantly  smaller  station  time  errors.  That  such  results  are 
achievable  when  locating  explosions  has  been  demonstrated  by 
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(among  others)  Chiburis  (1968)  and  Chiburis  and  Ahner  (1969)  and 
(197-0). 

However,  if  earthquakes  are  excluded^ from  analysis  of  location 
techniques  and  of  model  inhomogeneities,  the  number  of  regions 
available  for  study  is  severely  limited  and  the  results  are  apt  to 
be  less--general .  Also,  explosions  are  usually  detonated  in  such 
widely  separated  regions  that  the  large  variations  observed  in  the 
anomalies  between  regions  sometimes  cannot  be  assessed.  Since 
anomalies  at  a  given  station  are  known  to  change  significantly 
over  epicentral  distances  of  several  hundred  kilometers,  it  is 
almost  certainly  invalid  to  determine  anomalies  from  two  regions 
a  thousand  kilometers  or  more  apart,  linearly  interpolate  between 
them,  and  hence  obtain  anomalies  having  any  significance  for  under¬ 
standing  mantle  inhomogeneities. 

The  principal  difficulty  with  using  earthquakes  is  that  th^  - 
determination  of  ^  travel-time  anomaly  requires  an  event  location 
in  space:  latitude,  long itude ,. and  depth.  Knowledge  of  the  origin 
time  is  unnecessary.  Although  several  event  lists  are  available 
which  provide  hypocenters  from  which  anomalies  can  be  calculated 
for  any  network,  the  location  parameters  have  been  determined 
with  sets  of  stations  differing  in  number,  in  kind,  and  in  reli¬ 
ability  and  quality  o£  recording;  and  with  widely  different  earth 
models.  Thus  with  such  an  event  set,  the  station  anomalies  computed 
are  usually  not  consistent,  even  for  earthquakes  occurring  in  a 
single  region;  for  earthquakes  occurring  in  adjacent  regions,  the 
anomalies  agree  hardly  at  all.  For  explosions,  reproducibility 
has  been  verified  many  times  in  several  different  regions  (e.g., 
Chiburis  and  Ahner,  1970).  There  is  no  known  reason  (other  than 
mislocation)  why  earthquakes  should  yield  inconsistent  anomalies, 
but  explosions  consistent  ones. 

However,  a  study  of  earthquake  anomalies  need  not  be  deterred 
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if  one  is  aware  of  the  effects  that  mislocations  can  produce,  and 
the  conclusions  drawn  are  understood  in  that  context.  The 
objectives  of  the  present  study  are:  (1)  to  demonstrate  that 
location  and  anomaly  consistency  can  be  achieved  despite  nisloca- 
tion;  and  (2)  to  gain  improvement  in  the  relative  locations  of 
earthquakes  throughout  the  Aleutian  Islands  region.  Our  criteria 
for  success  are  location  and  anomaly  consistency. 

By  location  consistency,  we  mean  that  an  event  location 
obtained  with  a  well-distributed  network  must  agree  within  5  km 
with  the  locations  obtained  w'ith  well-distributed  subsets  of  the 

ri. 

initial  network.  It  must  also  be  true  that  the  least-squares 
time  residuals  of  the  several  network  solutions  for  the  event 
reduce  to  levels  generally  accepted  as  being  due  to  reading  and 
timing  error. 

* 

r» 

By  anomaly  consistency,  we  mean  that  the  anomalies  at  a 
station  must  change  across  the  Aleutian  Islands  in  a  physically 
plausible  way,  and  that  the  set  of  anomaly  functions  derived  for 
a  network  must  produce  location  consistency. 

The  sequence  of  steps  through  which  the  report  progresses  to 
achieve  the  desired  result  is  as  follows.  First,  it  is  demonstrated 
that  the  location  errors  obtained  by  using  different  networks  vary 
widely  in  magnitude  and  direction  when  locating  the  LONG  SHOT 
explosion  in  the  Aleutian  Islands. 

Therefore,  a  constant  network  is  selected  to  locate  a  set  of 
events  across  the  Aleutians  by  (1)  using  no  anomalies;  (2)  using 
the  LO.*G  SHOT  anomalies;  (3)  using  a  linear  anomaly  function  for 
each  station  (based  on  the  reported  locations);  (4)  using  a 
linear  anomaly  function  doubly  constrained  by  an  earthquake  in 
the  Fox  Islands  and  the  LONG  SHOT  explosion;  (5)  using  a  segmented 
linear  anomaly  function  constrained  by  four  events  equably  spaced  ® 


across  the  Aleutians;  and  (6)  repeating  the  last  step  bo  t  trans¬ 
lating  the  function  to  remove  the  observed  LO.\'G  SHOT  bias. 

Final 1> ,  instead  of  allowing  the  anomalies  to  varv  continuously 
between  constants,  the  anomaly  function  is  assigned  either  two 
or  three  constant  values  across  the  entire  Aleutians.  Each  of  the 
steps,  increasingly  more  complex,  provides  for  some  improvement 
in  location  consistency  and  all  are  presented  to  illustrate  the 
method  of  analysis  for  a  region  as  large  as  the  Aleutian  Islands. 

In  order  to  verify  the  results  obtained  above,  a  set  of 
synthetic  data  for  a  hypothetical  earth  model  is  analyzed;  it  is 
shown  that  the  same  order  of  consistency  is  achieved  when  the 
technique  is  applied  to  the  synthetic  data  as  it  was  to  the  real 
data. 


DESCRIPTION’  OF  THE  DATA 


\  set  of  we  1 1 - recorded  earthquakes  from  the  Aleutian 
Islands  region,  which  includes  the  Fox  Islands,  Andreanof 

t- 

Islands,  Rat  Islands,  and  Near  Islands,  was  selected  to  form 
our  data  base.  The  event  set  (Table  I)  is  composed  of  108 
earthquakes  and  two  explosions  ranging  in  magnitude  from 

a 

4.5  to  0.5.  The  earthquake  parameters  given  under  the  "NOS" 
column  are  those  reported  on  the  National  Ocean  Survey  (NOS) 

Preliminary  Determination  of  Epicenters  cards.  The  location 

\ 

parameters  of  LONG  SHOT,  an  underground  nuclear  explosion^  on 
Amchitka  Island  in  the  Rat  Islands,  are  those  released  by  the 
Atonic  Energy  Commission.  The  location  parameters  of  FLEXBAG, 
an  underwater  explosion  detonated  about  65  km  southwest  of 
Amchitka  Island  (Chiburis  and  Ahner,  1969,  for  a  seismic 
analysis) , 'are  those  reported  by  Kos  and  Kennedy  (1969).  A 
plot  of  thp  epicenters  is  shown  in  Figure  1.  / 

A  teleseismic  network  of  54  stations,  well  distribute4  in 
distance  and  azimuth,  was  selected  (Table  II);  of  these,  18  are 
World  Wide  Standard  Seismic  Stations  (WWSS  ),  28  are  Long  ftange 
Seismic  Measurement  (LRSM)  stations,  four  are  Geneva- type  Observa¬ 
tories,  and  four  are  F-ring  stations  at  the  Large  Ape^tuie  Seismic 
Array  (LASA)  in  Montana.  Arrival  times  were  read  for  those  events 
recorded  at  the  LRSM  stations  and  observatories;  arrival  times 
for  the  WWSS  stations  were  taken  from  the  Earthquake  Dita  Reports 
(EDR)  published  by  NOS. 
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I;i jjurc  1.  Ceographic  location  map  of  selected  Aleutian  Island  events 
Numbers  adjacent  to  locations  are  NOS  reported  depths. 
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DETER-*!  I  NATION  OF  TRAVEL-TIME  ANOMALY 


As  uscci  in  this  report,  the  travel-time  anomaly  at  station 
i  relative  to  station  j  is  defined  as 


where  I  is  observed  arrival  time  and  H  is  the  travel  tine  predicted 
iron  some  established  models.  This  definition  has  been  used  pre¬ 
viously  (e.g.,  Chiburis  and  Dean,  19^;  Chiburis,  I960,  1968).  In 
discussing  travel  tine  anomalies  for  a  given  network  of  stations,  we 
refer  travel  times  at  jill  stations  to  that  of  a  giveft  station 
of  ttie  network.  In  other  words,  i  varies  and  j  remains  constant. 

Travel-time  anomalies  were  computed  using  the  Herrin  (1968) 
tables  for  all  the  events  in  Table  I  and  for  those  stations  in 
Table  II  recording  the  events.  Initially,  the  input  latitudes, 
longitudes,  and  depths,  necessary  for  calculating  the  predicted 
times,  were  those  in  Table  I  reported  by  NOS. 

Plots  of  the  anomalies  for  two  selected  stations  relative  to 
UBO  are  shown  in  Figures  2a  and  2b  as  a  function  of  longitude 
across  the  Aleutians.  No  implication  is  made  as  to  the  longitude 
parameter  being  of  significance;  it  is  mere  convenience. 

The  variability  obseryed  in  Figures  ^a  and  2b  is  far  greater 
than  expected:  the  anomalies  are  known  to  change  as  a  function 
of  event  position,  but  not  in  the  erratic  manner  shown.  If  the 
results  in  these  two  Figures  were  representative  of  actual  anomaly 
'profiles,  the  problem  of  calibrating  sei smic' regions  to  the  degree 

necessary  for  accurate  location  work  would  be  almost  impossible. 

% 

If  it  is  assumed  that  the  station  anomalies  are  reasonably 
constant  within  a  region,  that  they  change  only  slowly  if  at  a(l 
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Figure  2a.  Anomalies  computed  from  NOS  hypocenters  for  station  NP-Ni. 
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fron  one  region  to  an  adjacent  one,  and  that  the  arrival  times 
are  read  correctly  to  within  0.2S  sec,  the  observed  erratic 
behavior  of  the  anomalies  can  be  attributed  only  to  mislocation 
of  the  hypocenter  from  which  travel-times  are  computed. 

It  is  well  established  that  the  anomalies  for  a  station  pair 
separated  by  a  hundred  kilometers  or  less  computed  from  events 
in  a  particular  region  are  constant,  and  that  the  anomalies  vary 
slowly  between  adjacent  regions  (Chiburis  and  Dean,  1965;  Chiburis 
1966,  196S ;  Chiburis  and  Ahner,  1969,  1970).  For  example,  in 
Figure  5  the  anomalies  at  stations  F4  and  E2,  relative  to  station 
AO  at  the  Large  Aperture  Seismic  Array  (LASA)  in  Montana,  are 
plotted  as  a  function  of  epicentral  distance  for  events  arriving 
along  a  northwesterly  azimuth.  In  this  figure,  the  anomalies  are 
virtually  constant  (within  about  0.1  sec)  for  a  particular  distance 
window  (anomaly  regiun)  and  vary  smoothly  from  one  region  to  the 
next.  The  effect  of  any  event  mislocation  in  this  example  would  be 
minimal  for  the  stations  used  because  of  the  small  distances  to 
the  reference  station  AO  (97  km  for  F4  and  6$  km  for  E2)  and  this 
fact  contributes  to  the  consistency. 

That  the  mislocation  effect  in  computing  the  anomaly  is  a 
function  of  station  pair  separation  has  been  discussed  by  Chiburis 
and  Dean  (1965)  and  by  Chiburis  1966.  In  the  latter  study  it  was 
shown  that  the  maximum  anomaly  error  6A  in  seconds  can  be  approximated 


a*6r 


where  a  is  the  distance  in  kilometers  separating  the  station  and 
its  reference,  5r  is  the  event  location  error  vector  in  km,  T  is 
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the  travel  time  in  seconds ,  and  A  is  the  epicentral  distance  in  km  to 
the  station  pair.  Therefore,  if  an  event  at  a  distance  of  5,000  km 
is  mislocated  by  50  km,  the  errors  in  the  computed  anomalies  for 
F4  and  E 2  would  only  be  about  0.10  sec  and  0.07  sec  respectively. 

On  the  other  hand,  if  the  station-pair  separation  is  4,000  km, 
as  is  the  case  for  XP-NT  and  UBO  in  Figure  2,  the  anomaly  error 
would  be  about  4  sec.  Clearly,  even  smaller  location  errors  (on 
the  order  of  10-20  km)  in  the  proper  direction  can  explain  quite 
large  disagreements  in  the  anomalies;  such  errors  are  the  probable 
causes  of  the  observed  anomaly  inconsistencies. 

Before  the  anomaly  errors  due  to  epicenter  mislocation  can  be 
properly  assessed,  the  depths  of  the:  events  must  be  known  as  well 
as  possible.  The  depths  resulting  from  an  unrestrained  least-squares 
solution  of  P  arrival  times,  without  corrections  for  travel-time 
anomalies,  are  known  to  have  an  average  error  of  75  km  or  more 
in  some  regions  (Chiburis  and  Ahner,  1970).  Such  errors  are  far 
too  large  to  permit  a  sensible  determination  of  travel-time  anomaly. 
The  only  known  way  to  determine  depth  unambiguously  is  to  identify 
the  phase  pP  and  restrain  the  least-squares  solution.  The  selection 
of  pP  on  an  event  sseismogram  recorded  at  a  single  station  is  highly 
subjective,  principally  because  of  signal-generated  noise  in  the  P 
coda  and  because  of  energy  from  a  multitude  of  other  phases. 

However,  if  a  suite  of  seismograms  recorded  at  different  stations 
which  are  widely  separated  in  distance  and  azimuth  is  aligned  ony 
the  P  wave,  time-distance  relationships  of  coherent  energy  can 
usually  be  correlated  with  distinct  phases.  Therefore,  110  events 
in  Table  I  were  analyzed  in  this  way.  The  maximum  disagreement 
with  the  NOS  Reported  depths  in  the  event  set  is  about  50  km. 
Although  depth  errors  by  themselves  do  not  produce  large  anomaly 
changes,  the  differences  in  epicenters,  had  the  events  been 
located  with  restrained  depths,  can  .be  large  and  can  produce 
significant  anomaly  changes.  Therefore,  all  of  the  events  in 
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Table  I  were  relocated  with  depths  restrained  to  our  values  based 
on  pP,  but  Using  the  same  stations  as  were  used  by  NOS  and  reported 
in  the  EDR.  The  number  of  stations  varies  between  12  and  126, 
excluding  LONG  SHOT.  The  resulting  restrained  epicenters  are  given 
under  the  ”Adj us  ted”  column  in  Table  I.  Also  included  are  columns 
indicating  the  estimated  reliability  of  the  pP  time  pick  and  the 
shift  vectors  from  the  NOS  locations  to  the  Adjusted  locations. 

The  maximum  shift  is  about  40  km,  excluding  the  nine  events 
appearing  in  Table  Ic.  (These  nine  events  are  unusual  in  their 
signal  characteristics  and  in  their  anomalies;  it  is  believed  that 
they  are  multiple  events,  and  they  will  be  analyzed  in  a  special 
study  at  a  later  date.)  The  reasons  for  the  fact  that  there  are 
slight  differences  between  the  number  of  stations  used  by  SDL 
and  by  NOS  are  as  follows:  (1)  no  stations  were  included  which 
reported  PKP  arrival  times  or  two  P  wave  arrivals  close  together; 

(2)  where  two  stations  are  located  at  virtually  the  same  site 
(e.g.,  MBC  and  NP-NT)  only  one  was  used. 

Recomputing  the  anomalies  from  the  Adjusted  locations  for 
stations  NP-NT  and  NUR  yields  the  anomaly  plots  shown  in  Figures 
4a  and  4b.  A  comparison  of  these  results  with  those  in  Figures  2a 
and  2b  shows  essentially  no  improvem'  t  in  the  functional  patterns. 

If  it  is  assumed  that  some  of  the  scatter  is  due  to  those  events 
not  lying  within  thq  Aleutian  Islands  region  proper,  these  events 
can  be  separated  and'  so  noted  by  the  circled  values  in  Figures  4a 
and  4b  .  A  total  of  2J.  events  was  eliminated  in  this  way;  they  are 
listed  in  Table  lb.  However,  the  anomalies  still  do  not  show  a 
consistent  enough  rel 
purposes. 

The  question  of  reading  errors  being  the  producer  of  anomaly 
errors  can  be  only  partly  answered.  Generally,  for  the  LRSM  and 
V£LA  stations  listed  in  Table  II  and  the  examples  in  Figures  2  and  4, 


tion  across  the  region  to  be  used  for  location 
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.figure  4b.  Anomalies  computed  from  adjusted  hypoccnters  for  station 


all  arrival  times  were  read  by  SDL  analysts^/so  reading  errors 
are  believed  to  be  acceptably  low.  Arrival ''times  for  other 
stations  were  taken  from  the  Earthquake  Data  Reports  published 
by  the  NOS.  Later  in  this  report  we  show  that  these  reported 

times  are  often  in  error  by  as  much  as  one  seconds 

\ 

-  j 

Therefore,  the  event  set  in  Table  I,  although  largely 
composed  of  well-recorded  earthquakes,  most  probably  contains 
location  errors  of  perhaps  30  km,  or  more  and  the  anomaly  incon¬ 
sistency  is  not  surprising.  The  possible  errors  one  might  encounter 

•  *  / 

when  locati/ng  events  in  the  Aleutian  Islands  is  demonstrated  in 
/  ~ 
the  next  section. 


ERRORS  IN  LONG  SHOT  LOCATION  USING  DIFFERENT  STATION  NETWORKS 

<  '  " 

The  underground  nuclear  explosion  LONG  SHOT  was  detonated  * 

29  October  1965  on  Amchitka  Inland  in  the  Rat  Islands  with  a 
yield  of  80,000  tons  TNT  equivalent.  P  wave  signals  Were  well 
recorded  at  virtually  all  distances  and  azimuths.  In  their  com-  ~ 
prehensive  report,  Lambert  e^  al.,  (1969)  discuss  the  seismic 
location  made  with  travel-time  data  from  329  stations  available 
to  them,  although  probably  several  hundred  more  stations  recorded 
the  event.  The  location  obtained  by  restraining  the  focal  depth 
to  the  known  value  and  using  either  the  Herrin  (1968)  or  Jeffreys  - 
Bullen  tables  was  approximately  21  km  in  error  to  the  northwest* 
(Lambert  et  al.,  1969).  To  now  demonstrate  the1  possible  variations 
in  location  for  the  Rat  Islands  region,  selected  subsets  of  the 

329-station  network  are  taken  in  different  ways  and  the  resultant .. 
locations  discussed.  , 

First,  a  network  can  be  defined  on  the  basis  of  its  .izimuth 
aperture  or  quadrantal  coverage.  Requiring  a  minimum  aperture  of 
180°,  or  two  quadrants,  and  beginning  with  the  sector  0°-180°  (north¬ 
east  and  southeast  quadrants),  all  stations  within  that  azimuth 
range  regardless  of  distance  are  used'to  locate.  Incremr  ■ ing 
the  two-quadrant  sector  by  22.5%  a  total  of  16  separate  locations 
can  be  made  with  networks  all  having  180°  aperture.  These  results 
are  shown  in  Figure  5,  where  the  329-station  case  is~included  for 
comparison,  as  well  as  the  rmmber  of  stations  for  each  case,  which 
varies  from  57  to  27>STrafions .  The  resultant  location  errors  are 
as  small  as  13  km  and  as  large  as  45  km.  There  appears  to  be  a 
preferred  shift  toward  the  northwest,  suggesting  a  bias  independent 
of  the  network  although  five  of  the  16  locations  shift  otherwise. 

From  the  results  of  this  Figure,  one  can  state  with  some 
assurance  that  the  epicenters  of  events  in  the  Rat  Islands  and 
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Andreanof  Islands  (and  probably  throughout  the  entire  Aleutian 
Islands)  are  not  known  to  within  20  km  of  their  true  locations, 
regardless  of  the  distribution  or  stability  of  the  networks  used 
to  locate  them.  A  reiteration  of  a  conclusion  previously  drawn 
(Chiburis  and  Ahner,  1970)  is  in  order:  location  errors  are 
caused  by  the  set  of  travel-tine  anonalies  at  those  stations 
uniquely  defining  a  particular  network;  several  networks,  similarly 
distributed  but  composed  of  different  stations,  can  yield  entirely 
different  locations  for  the  same  event,  even  if  a  randomly  selected 
network  sometimes  shows  a  bias  due  to  the  source. 

If  single-quadrant  distribution  (90°  aperture)  is  accepted 
as  a  minimum,  instead  of  two  quadrants ,  and  the  quadfant  is  incre¬ 
mented  by  22.5°, the  resultant  locations  are  shown  in  Figure  6^ 

The  errors  now  range  between  2  km  and  163  km,  with  an  average  of 
more  than  50  km  and  no  apparent  preferred  shift  direction.  Tb'e 
number  of  stations  varies  between  22  and  200  with  no  firm  corre¬ 
lation  with  the  size  of  the  error. 

If  the  entire  network  of  329  stations  is  now  ordered  on  the 
basis  of  increasing  azimuth  (0°  to  360°),  and  every  tenth  station 
of  the  ordered  set  is  used  for  locating,  a  set  of  nine  32-station 
networks  is  obtained,  each  with  approximately  four-quadrant 
coverage .  The  location  results  are  shown  in  Figure  7,  where  a 
distinct  northwesterly  error  of  about  17  km  can  be  observed. 

This  suggests  that  for  the  earth  model  used  to  compute  travel 
times,  those  stations  that  "bear  northerly"  from  LONG  SHOT  are 
early  in  time  relative  to  those  stations  that  "bear  southerly". 

A  plot  of  the  zero-mean  travel-time  anomalies  vs  azimuth  is 
shown  in  Figure  8,  where  the  northerly  stations,  from  340°  to  70° 
azimuth,  are  indeed  seen  to  be  generally  early  and  the  southerly 
stations  from  70°  to  340°  generally  late,  However,  the  scatter  of 
the  anomalies,  even  within  azimuth  windows  of  20°-30°, approaches 
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Figure  7.  LONG  SHOT  location  errors  using 
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increasing  azinuth. 
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3  sec,  which  in  itself  is  enough  to  cause  location  errors  of 
50  kra.  Because  of  this  fact,  the  practice  of  fitting  a  sinusoidal 
curve  through  these  data  in  order  to  predict  what  the  anomalies 
would  be  for  events  occurring  in  this  region  is  hardly  more  than 
an  exercise,  although  It  may . remove  a  portion  of  the  bias.  For 
precise  location  work,  anomalies  need  to  be/  known  far  better  than 
0.5  sec;  but  with  a  known  scatter  of  3  sec,  any  reduction  in 
location  errors  obtained  by  applying  sinusoidal  azimuth-dependent 
corrections  is  entirely  coincidental.  If  instead  of  ordering  the 
station  set  by  azimuth,  ordering  is  on  the  basis  of  increasing 
epicentral  distance  and  every  10th  station  selected,  another  set 
of  nine  32-station  networks  is  obtained,  again  each  with  approx¬ 
imately  four-quadrant  coverage.  The  results  are  shown  in  Figure  9, 
which  show  about  the  same  northwesterly  bias  as  in  the  azimuth- 
ordered  results  of  Figure  7. 

In  any  least-squares  procedure  of  location,  the  distribution 
of  the  true  (zero-mean)  time  errors  regardless  of  distance  or 
azimuth  is  critical  to  the  result.  In  order  to  solve  for  the 
correction  coefficients,  the  residual  errors  are  assumed  to  have 
the  normal  distribution  about  an  assumed  travel-time  distance 
relationship.  But  the  actual  errors  for  a  particular  event,  nay 
not,  and  indeed  generally  do  not,  have  the  normal  distribution. 

A  plot  of  the  329  time  errors  computed  from  the  true  location  of 
LONG  SHQJ  is  shown  in  Figure  10a  for  a  1.0  sec  error  interval. 

There  is  seen  to  be  a  definite  skewness  toward  negative  errors  or 
early  arrival  times  relative  to  the  earth  model;  the  standard  devia¬ 
tion  for  these  data  is  1,3  sec.  This  is  another  effect  of tShe  varia¬ 
tion  of  travel  time  with  azimuth.  After  least-squares  adjusting  the 
epicenter,  resulting  in  a  21  km  shift,  the  minimized  zero-mean  time 
errors  are  shown  in  Figure  10b.  Now  the  nearly-normal  distribution 
is  clear,  the  standard  deviation  for  these  data  being  1.0  sec.  The 
assumption  made  by  employing  a  least -squares  scheme  in  the  first  place 
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Figure  9.  LONG  SHOT  location  errors  using  four-quadrant,  32-station 
networks  conposed  of  every  tenth  station  selected  by 
increasing  distance. 
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is  that  any  subset  of  travel-time  anomalies,  even  a  large  one  (Figure 

10a),  has  a  normal  distribution  and  must  be  minimized  because  of 

/ 

random  imperfections  not  because  of  systematic  early  or  late 
arrival  times.  The  anomalies  are  due  to  lateral  and  vertical 
inhomogeneities  which  are  not  allowed  for  in  the  earth  model; 
the  actual  effects  of  t,hese  inhomogeneities  may  be  large  and 
nothing  dictates  that  a  particular  set  of  anomalies  for  the 
network  locating  the  event  should  be  as  small  as  possible  in  a 
least-squartes  sense.  The  only  errors\  which  should  ever  be 
minimized  are  those  due  to  reading  precision;  these  are  random 

and  probably  have  a  zero-mean  normal  distribution. 

'•  \ 

From  the  above  results,  an  example  can  be  given  of  the 
difficulty  in  working  with  a  mislocated  earthquake  data  set  and 
of  possible  anomaly  errors  computed  from  the  various  LONG  SHOT 
solutions , for  station  XP-XT.  The  anomaly  inconsistency  for  this 
station  is  shown  in  Table  III,  in  which  the  range  is  11.44  sec. 

If  LOXG  SHOT  had  been  an  earthquake,  any  one  of  these  epicenters 
might  have  been  reported,  depending  on  the  particular  network 
available  at  the  time  of  occurrence,  and  hence,  any  one  of  these 
anomalies  would  have  been  computed  aind  accepted.  Out  if  accurate 
locations  are  to  be  realized  for  a  particular  region,  anomalies 
must  be  known  to  within  0.25  sec,  including  reading  error;  an 
11  sec  uncertainty  is  useless. 

The  principal  conclusion  to  be  drawn  from  the  preceding  ■  , 

exercise  is  that  by  studying  anomalies  across  a  region  as  large 
as  the  Aleutian  Islands  (approximately  1600  km  in  extent)  and 
that  by  using  earthquakes  which  are  all  differently  mislocated, 
the  anomalies  must  be  expected  to  display  variations  and  inconsis¬ 
tencies  which  are  not  real^  However,  even  if  event  mislocation 
is  the  chief  cause  of  the  Variations,  an  attempt  can  be  made  to 
achieve  both  anomaly  and  location  consistency,  no  matter  how 
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artificial  the  anomaly  functions  and  the  relative  locations 
appear  to  be.  With  the  LONG  SHOT  explosion  providing  a  -bias 
constraint  in  the  vicinity  of  Amchitka  Island,  relative  loca 
tions  can  be  made  using  functional  anomalies  at  increasing 
distances  from  LONG  SHOT,  always  demanding  least -squares 
consistency  and  reasonable  functions.-. 
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DETER?! INAT ION  OF  ANOMALY  FUNCTIONS 

Selection  of  a  constant  network 

It  uas  shown  in  a  previous  report  (Chiburis  and  Ahner,  1970)*, 
as  well  as  in  the  preceding  section,  that  "location  bias"  is  a 
function  of  the  particular  network  used;  also,  if  a  constant 
network  is  used  for  locating  a  set  of  events  in  one  region, 
relative  accuracv  within  the  event  set  is  not  affected,  although 
all  of  the  events  nay  be  translated  by  an  unknown  bias.  However, 
if  different  networks  are  used  to  locate  different  events,  the 
results  of  locating  in  a  particular  region  display  large  unpre-. 
dictable  errors,  unless  anc  tlies  are  applied. 

Therefore,  in  this  report,  in  order  to  remove  the  network 
effect,  a  constant  network  is  selected,  all  stations  of  which 
recorded  a  suitable  number,  of  the  events,  in  T^ble  I.  Subsets  of 

\  J  •  /'  •  • 

this  network  caH  then ’be  used  to  demonstrate  location  consistency. 

Initially,  the  constant  network  was  composed  of  seven  LRSM 
and. VELA  stations  and  eight  NOS  and  participating  stations.  The 
number  of  events  was  ten.  Various  techniques  were  initially 
employed  to  achieve  consistency, .none  of  which  yielded  satis¬ 
factory  results.  \Eventt^a^l>*  it  was  supsected  that  the  reported 
arrival  times  at  ?vOS  stations  were  partly  responsible  for  the 
poor  results.  An  effort  was  made  to  obtain  the  seismograms  from 
the  NOS  library  in  order  to  check  the  readings.  Only  those  seis¬ 
mograms  from  the  KKSS  system  were  available;  the  seismograms  from 
participating  stations  we^^oot.  Table  IV  is  a  list  of  only 
those  reading  discrepancies  greater  than  1.0  sj»c  found  at  eleven 

.  U  T. 

selected  KfcSS  stations  for  a  set  of  50  events. 

.  ''j 

c  A  total  of  2“9  readings  was  compared.  Of  these,  40  were  different 
by  Adre  than  1.0  sec  and  114  by  more  than  0.5  sec.  Consequently, 
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TABLE  IV 

Selected  Arrival-Time  Discrepancies  Greater  than  1.0  Sec 

(SDL  -  NOS) 


-1.4  +  1.8  -1.7  -1.1  -2.0  -1.0  -2.0  ♦1.4 

+1.0  +1.4  +1.4  -1.5  -1.9  +1.7  +2.1  +1.4 

+2.1  +1.0  -1.0  -1.2  -1.0  -1.0  -1.5  -1.8 

-1.4  -1.7-  -2.4  -1.1  -1.0  +1.1  -2.1  -1.7 

-1.0  +1.0* -1.0* -1.1  -1.3  -19.9*  +1.2  -1.3 


♦Reading  not  used  by  NOS  in  location  reported, 
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an  readings  used  in  the  rest  of  this  report  are  only  those  I 
for  which  seismograms  were  obtained  and  analyzed  at  SDL. 

The  constant  network  finally  selected,  referred  po  as 
Network  1.  is  composed  of  fourteen  stations  Three  ^sets 
this  network  were  devised  (Table  V).  Also 

are  seventeen  events,  each  of  which  was  recorded  ^  all  four 
teen  stations.  The  LONG  SHOT  explosion  is  included 

but  „ot  FLEXBAG,  because  it  was  not  recorded  by  all  four 

^  stations .  Therefore,  only  one  event  serves  as  bias  con  ro 

,  .  r  seventeen  events  is  snown 

(nr-  the  study.  A  geographic  map  of  the 

:F;  :r:  IK  Although  farther  apart  than  desired,  the  number 

Of  events  is  the  maximum  obtainable  for  as  many  as  f“urt®en  . 

stations  By  decreasing  the  number  of  stations,  one  could  obtain 

«  but  the  network  would  be  less  effective  in  azimuth 

and  distance  distribution. 


!n; Lial  locations  \ 

Using  the  Adjusted  locations  given  in  Table  I  as  the  input 

r  the  seventeen  events  were  located  by  the  four  networks 
Pa"' ””  a;rlval  times.  Recall  that  the  Adjusted  locations  are 

relocations  of  NOS  hypocenters  with  depths  restrained  to  pR  readings 

» -“irr  r= 

"!H  i....  ' 

scatter  within  the  clusters  is  as  large  as  different 

for  this  poor  result  is  that  the  networks  are  composed  of  different 

sets  of  stations  which  have  quite  different  distributions 

i •  The  clustering  and  standard  deviations  of  the  leas 
anomalies.  The  clustering  listed  in  Table  VI. 

squares  solutions  for  the  seventeen  events  are  listed  in  T 

Ceveral  standard  deviations  become  quite  small,  showing 
^1  resides  are  a  necessary  but  not  ,  suffici.n,  indication  of 


accuracy. 
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Figure  11.  Geographic  location  map  of  17  selocted  evonts. 


LONG  SHOT  anomalies  for, corrections 


The  LONG  SHOT  anomalies,  which  are  accurately  known,  are  given 
in  Table  VII  for  the  fourteen  stations.  If  it  is  assumed  that  these 
anomalies  are  valid  across  the  entire  Aleutian  region,  they  can 
be  applied  as  corrections  to  all  arrival  tines.  The  location 
results  of  doing  this  are  shown  in  Figure  13.  Several  of  the 
events,  particularly  in  the  Fox  Islands,  now  tend  to  approach 
a  tighter  cluster  of  solutions  for  the  four  networks.  But  for 
the  Andreanof  9,  10,  11,  and  the  Rat  Near  events,  no  improved 
clustering  is  observed.  In  fact,  when  no  anomalies  were  applied, 
far  greater  clustering  was  obtained  for  these  events  (Figure  12); 
regardless  of  the  clustering,  the  standard  deviations  are 
generally  now  larger  (Table  VIII).  This  implies  that  a  poorer 
travel-time  fit  is  being  obtained  in  the  least-squares  scheme 
when  LONG  SHOT  anomalies  are  used  as  corrections. 

Both  network  clustering  and  a  reduction  in  the  standard 
deviation  (or  an  approach  to  an  acceptable  value  --  say  0.25  - 
0.50  sec)  are  necessary  in  order  to  achieve  consistency  and 
relative  accuracy;  but  they  are  not,  alone  or  together,  sufficient 
to  obtain  absolute  accuracy  (complete  removal  of  "bias").  The 
only  known  way  to  achieve  absolute  accuracy  is  to  have  seismic 
data  available  from  a  known  explosion  previously  detonated  in 
the  vicinity.  One  might  expect  that  LONG  SHOT  would  approximately 
calibrate  adjacent  regions.  In  fact  toward  the  east  the  standard 
deviations  decrease  for  Andeanof  -11,  -10,  and  -9  when  using 
the  LONG  SHOT  anomalies.  However,  the  clustering  results  for 
these  same  events  are  worse;  thus  the  results  are  inconclusive. 

As  the  LONG  SHOT  anomalies  are  not  valid  except  in  the 
immediate  vicinity  of  Amchitka  Island,  and  possibly  eastward  for 
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TABLE  VII 

LONG  SHOT  Anoaalies  for  14  Stations  Relative  to  UBO 


Station 

NDI 
QUE 
NUR 
STU 
NOR 
XP -XT 
COL 
HX-ME 
RK-OX 
LAO  10 
NMO06 
UBO10 
TF060 

rue 


Observed  Long  Shot 
Anonalies  in  Seconds 

-1.23 

-0.10 

-1.4(> 

-2.74 

-0.72 

0.02 

-1.39 

-2.21 

-2.12 

-1.29 

-0.S4 

0.00 

0.62 

0.10 


0 


1 


Location  shifts  for  four  networks  obtained  by 
applying  LONG  SHOT  anomalies. 


a  short  distance,  a  set  of  anomaly  functions,  one  for  each  station, 
needs  to  be derived  which  will  provide  consistency ,  ifytwork 
clustering,  and  acceptable  standard  deviations  aerq^s  the  entire 
Aleutian  Islands.  It  would  be  trivial  to  draw  functions  exactly 
through  the  anomalies  in  Figures  4a  and  4b,‘  which  fire  simply  the 
anomalies  from  the  Adjusted  locations,  for  two  reasons:  first, 
it  is  hoped  that  the  anomalies  would  not  change  as  drastically 
as  shown  over  such  small  distances,  and  second,  that  there  should 
be  at  least  a  suggestion  of  the  anomaly  being  monotonic  or  having 
a  smooth  rate  of  change  over  distances  of  500  km  or  so.  In  other 
words,  the  functions  tha^  one  derives  should  be  geophysically 

realistic.  ' 

\ 

Simple  functional  anomalies  , 

If  it  is  assumed  that  the  anomaly  variations  produced  by  the 
earth  are  greater  (than  the  errors  due  to  dislocation  or  to 
measurement,  the  Adjusted  locations  can  be  used  to  compute 
anomalies  at  all  of  the  stations  and  to  estimate  anomaly  functions 
by  fitting  a  linear  relation  through  the  anomalies  for  each  station. 
Because  the  LONG  SHOT  anomalies  are  accurately  known,  bias  con¬ 
straints  are  available  for  each  of  the  curves  on  the  anomaly  axes.. 
These  curves  and  their  translations  are  shown  in  Figures  14a  and 
14b. 

A  linear  function  is  the  simplest  possible  and  provides  a 
starting  point  from  which  more  complex  functions  can  be  devised-. 

If  these  values  are'now  applied  as  a  function  of  longitude,  the 
event  set  can  be  located  by  the  four  networks,  the  results  of  which 
are  plotted  in  Figure  15  and  listed  in  fable  IX.  Except  for 
Andreanof- 11,  it  cah.be  seen  that  both  clustering  and  reduction  in 
standard  deviation  ate  being  accomplished  about  the  same  as  in 
Figure  13.  The  four  network  solutions  for  an  event  aie  now 
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t  •  i • .  . .1 :  :  ieviat  ions  are  remarkably  improved  through  the  Fox 
\:.d:  t:.c:  islands,  generally  hy  a  factor  of  about  two.  There- 

,  ,:s :  ;.g  two  widelv  separated  events  as  constraints,  the  least- 
s,..ares  fits  and  the  network  clustering  of  events  between  the  two 
i:e  ette:  than  if  no  anomalies,  constant  anomalies,  or  single- 
.  n  t:  ii.nt  linear  functional  anonalies  had  been  applied. 

in  order  to  gam  location  control  in  the  Near  Islands,  and 
to  improve  the  locations  in  the  Andreanof  Islands,  two  additional 
c  astrair. :>  can  be  imposed,  one  in  each  of  tlie  two  regions. 

•  AR.\  anomalies  . 

neca use  the  anomalies  are  neither  constant  nor  linear  across 
tile  entire  Meutian  region,  a  set  of  anomaly  functions  can  be 
derived  winch  are  (Ij  constrained  by  the  anonalies  computed  from 
four  assumed-correct  events  (Fox-1,  Andreanof-6,  LONG  SHOT,  and 
'.ear-1")  and  id  piecewise  linear  between  the  four  constraints. 

A  plot  of  the  functions,  referred  to  as  "FARM”  (for  £ox-Andreanof - 
Rat -Near j,  is  shown  in  Figures  ISa  and  ISb. 

Ihe  location  shifts  obtained  by  applying  the  FARN  anonalies 
are  shown  in  Figure  19.  The  scatter  of  the  clustered  locations 
for  each  event  is  reduced  and  is  now  about  the  sane  across  the 
entire  Aleutian  region,  except  for  Near-15  and  Andreanof- 11 .  Even 
for  these  two  events,  if  Network  4  is  not  included,  the  clustering 
approaches  the  average.  It  can  be  expected  that  Network  4  will  be 
unstable  due  to  its  snail  aperture  of  about  59°  where  the  effect 
of  snail  reading  errors  is  exagj^rated .  A  sunnarv  of  the  results 
is  given  in  Table  .X4p~i?he re  all  of  the  standard  deviations  becone 
less  than  0.5  seconds,  nost  being  less  than  0.35  seconds.  The 
clustering  error  is  about  8  kn,  and  without  Network  4  being 
considered  for  Andreanof-11  and  Near-15,  it  is  about  6  kn. 


Figure  13a.  Deterninat ion  of  anomalies  using  Fox,  Andresnc 
and  Near  fF.ARN)  constraining  events. 
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Rat  Island  event,  the  anomalies  from  this  fiducial  location 
in  error;  can  be  computed  and  compared  tc  the  anomalies  from 
Adjusted  seismic  locations  c i  earthouakes.  -i  ^et  o* 
functional  anomalies  can  then  be  fitted,  without  con 
all  of  the  observed  anomalies.  The  or.lv  reoui  remer.t  cm  the  * 
is  that  they  be  piecewise  linear;  the  n 

derived  for  all  of  the  stations  varies  from  two  to  three, 
anomaly  plots  and  the  fiducial  functions  for  the  station; 
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:r  ;_;„re5  22a  and  22b.  The  observed  true  LONG  Sli  i  anonalies 
j.-'-j  included  for  comparison.  Significantly,  for  every  station,  tie 

. _ ;.  si;  7  anonalv  agrees  much  better  than  the  true  ar.otai.. 

„  • ti.e  earthquake  anomalies.,  which  supports  the  hypothesis  that 
- v  Aleutian  Islands  earthquakes  in  the  vicinity  of  LONG  SHOT  are 
_:sl;c.sted  ;r.  a  similar  way. 

f  the  functions  in  Figure  20  are  now  applied  as  corrections, 
t  .e  locations  obtained  with  the  fou.  networks  are  shown  in  Figure  21 
^r.d  the  results  summarized  in  Table  XII. 

For  nine  of  the  events.  Network  4  yields  a  result  clearly  at 
..st  ance  with  the  other  networks;  for  Networks  2  and  5,  there  is 
one  event  at  variance.  Therefore,  the  results  for  these  cases 
acre  .net  included  in  forming  the  average  shift  from  the  mean  or 
:r.  fcrmir.c  the  "best  estimate"  of  the  epicenter.  The  excluded 
results  are  marked  with  an  asterisk  on  Table  XII. 

jye  average  shift  from  either  the  ne an  or  "best  estimate  is 
about  5 .  S  km,  and  for  a]_l  networks  the  standard  deviations  are 
less  than  0.5  sec  (except  for  one  result,  the  standard  de\iations 
arc  all  less  than  0.54  sec J  with  the  average  being  less  than  0.2a 

an  acceptable  value  to  indicate  reading  precision.  He  should 
.note,  hcwcF>s-r  that  if  the  events  with  asterisks  are  included,  then 
the  deviation  from  the  mean  is  S.5  kilometers,  about  the  sane  as 
for  the  FASN  anomalies. 

if  anomalies  are  now  computed  from  the  mean  epicenters 
derived  above,  the  results  as  shown  in  Figures  22a  and  22b 
indicate  that  the  station  anonalies  display  less  scatter  fron 
a  function  than  before  the  shifts. 

It  is  now  necessary  to  correct  the  anomaly  functions  for  the 
known  LONG  SHOT  bias  and  to  attempt  to  smooth  them  (remove  sharp 
changes )  across  the  Aleutian. 
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|;iKUrc  21.  Location  shifts  for  four  networks  obtained  !>' > 

fiducial  l.ONd  SHOT  functional  anomalies* — 


Mcutian  functional  anonalies 

Because  it  is  known  that  LONG  SHOT  is  actually  r.iblcta.ei 
bv  21  kn  to  the  north-northwest,  we  see  that  to  obtain  anonalies 
wnich  yield  snail  absd.ute  location  errors  LONG  SHOT,  and  to 
sone  degree  the  other  events,  must  be  shifted  toward  the  .south¬ 
east.  f  ne  would  choose  to  shift  then  in  such  a  way  that  the 
ario.naly  functions  are  further  sncothed  out.  One  cl  tne  guiding 
observations  in  determining  which  events  to  sr.il  t,  and  oy  how 
nuch,  is  the  apparent  "anomaly  discontinuity"  displayed  by 
nany  of  the  stations  (especially  STU,  NDI.  ,’’n,  and  Rh-  NJ  in 
the  vie  ini  tv  of  179'K  (Andreanol -10  and  -llj.  ihe  rcaso^  -cr 
this  discontinuity  is  traceable  to  the  original  nislcoation 
effect  in  the  Aleutians,  but  why  it  exists  is  presently  not  known. 
Anotner  guiding  observation  is  that  the  only  quantitative  bias 
information  available  is  that  provided  from  the  LONG  SI.  I  shift, 
which  amounts  ro  21  kn  toward  the  north-northwest.  Kith  these 
guidelines,  we  choose  to  shift  the  positions  of  all  of  the  events 
between  .Andreanof - 1 1  and  Near-1"  21  kn  to  the  south- southeast 
from  their  respective  "best  estimates".  All  events  between 
Andreanof- 10  and  Fox-1  were  not  shifted  at  all  from  their  "best 
estimates".  The  plots  of  the  anonalies  computed  from  the  variably 
shifted  locations  and  the  new  anomaly  junctions  which-re&ult  are 
shown  "in  Figures  23a  and  25b-  Although  the  discontinuity  remains,, 
it  has  been  reduced,  and  the  functions  are  slightly  more  continuous 
Ke  choose  these  "Aleutian  functional"  anomalies  to  be  the  best 
ones  to  apply  when  any  combination  of  stations  in  Network  1  is 
used  to  locate  events  in  the  Aleutian  Islands  region. 

o, 

The  results  obtained  by  applying  the  Aleutian  functional 
anomalies  are  shown  in  Figure  24  for  the  four  networks  and  are 
summarised  in  Table  XIII.  The  corresponding  final  "best  estimates" 
of  the  epicenters  are  given  in  Table  XIV  with  the  shifts  from 
both  the  original  NOS  and  the  Adjusted  locations, in  .able  I. 
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Determination  of  Aleutian  functional  anonalies 
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j:  jii  ure  24 „  Location  shifts  for  four  networks  usinK  Aleutian 
functional  anomalies. 


TABLE  XIV 

Final  Best-Estinates  of  Event  Set  Located 
vith  Aleutian  Functional  Anoaalies 


Fvent 

Mean  or  3est-Estinate 
Locations 

Latitude  Longitude 

fnerreesl  (Degrees)  _ 

Shift 
XOS  LOC 

to 

Fox 

1 

53.583  X 

165.727  K 

12 

Fox 

52.-43  X 

166.931  V 

^  1 

Fox 

3 

52.646  X 

168.163  K 

6 

Fox 

4 

52.082  X 

171.372  K 

8 

Fox 

3 

'  52.008  X 

1-1.535  Sf 

19 

And 

* 

51.682  X 

1-73.361  V 

8 

And 

i 

51.609  X 

173.487  K 

11 

And 

a 

51.402  X 

1^6.121  V 

10 

And 

9 

51.452  X 

178.482  K 

\  J 

5 

And 

10 

51.528  X 

1-8.385  K 

4 

And 

11 

51.088  X 

179.628  X 

30 

Long  Shot 
Rat  13 
Rat  14 
Near  15 
Near  16 
Near  17 


51.472  X 
51.703  X 
51.663  X 
52.292  X 
52.670  X 
53.004  X 


179.214  E 
176.507  E 
176.196  E 
1731621  E 
172.114  E 
171.217  E 

Average  Shift 


5 

14 

28 

17 

28 

12 

13 


Shift  From 
Adjusted  Location 
_ to _ 

8 

5 

10 

3 

6 
1 

1  8 
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3 

10 

26 

5 
9 

27 

19 

24 

17 
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Several  things  have  been  accomplished  with  these  Aleutian  ^ 

..  rn  a  tisht  location  cluster  (within  aJOJt 

functional  anomalies,  (lj  a  tign 

iin  8  In  if  one  includes  the  readings  with  astens  s  us  «• 
three  or  four  different  networks;  (2)  acceptably  low  standard 
deviations  (aoout  0.25  sec)  for  all  networks  and  events;  an.  ), . 
reasonable  functional  anomalies  for  eacn  station. 

The  question  persists,  though,  of  the  absolute  ar.c  rela‘“c 
accuracv  of  location  throughout  the  entire  region.  ^  ne  anor.au, 
functions  were  derived  using  assumed  locations  for  tr.e  even  s._ 

That  these  assumed  locations  are  incorrect  by  an  unk..o..„  a .  ou 
is  beyond  reasonable  doubt.  However,  as  'ong  as  it  is  w.  .a  -■  ^ 

events  were  used  to  determine  the  functions,  all  5uD>cqjem  e.e..t. 
can  be  referenced  to  them,  almost  as  if  it  did  not  matter  -a 
their  actual  location  errors  may  have  been,  alien  and  if  an  event 
occurs  somewhere  in  the  Aleutian  Islands  other  than  on  Anchi  ha 
island  under  conditions  such  that  its  location  is  accurately 
known  ’the  anomalies  can  be  recomputed  and  the  anomaly  tunc t tons 
translated  to  accommodate  the  additional  constraint.  ine  posi.t  - 
of  this  additional  event  may  be  accurately  known  due  to  it 
another  explosion  (far  away  from  Anchitka) ,  due  to  its  being 
recorded  bv  a  local  high  quality  network,  or  due  to  clear  sur  .. 
intensity  effects  (faulting,  jointing,  etc).  Th.  latter  criterion 

can  only  serve  to  remove  gross  bias  ^ distance 

’J( from  earthquakes  depend  on  man?  factors  P 

this  is  further  complicated  by  the  Aleutian  Islands  region  itself, 
being  composed  of  small  and  sometimes  uninhabited  is.ar.ds  ,..  tnereo.. 
limiting  the  areal  observation  of  intensity,  out  regaruless  . 
the  availability  of  calibration  events,  as  pointed  out  ,n  an 
earlier  section  it  is  important  that  one  first  achieve  cons  te..  . 
throughout  a  region,  and  then  concentrate  on  removing  bia.  b. 
detailed  enpi  ’leal  studies. 


liic  nature  of  these  studies  can  be  briefly  outline^.  Sclent 
a  network  composed  of  nore  stations  than  were  used  in  this  report, 
perhaps  50,  such  that  at  least  six  or  eight  stable  network  subsets 
can  be  used.  Unstable  networks  (as  Network  4)  should  not  be  included 
increase  the  number  of  events  to  a  far  larger  nunber  tnan  used  here, 
perhaps  40.  Critically  read  and  re-read  arrival  tines  for  all 
events  at  all  stations.  Using  the  constraints  of  network,  clustering, 
goodness -of - f i t ,  anomalies  fron  known  explosions  or  "well  located 
earthquakes,  reasonable  functions  (perhaps  linear,  but  at  least^ 
smoothly  carving),  and  naxinura  permissible  anomaly  values,  one 
could  program  and  determine  in  a  least-squares  sense  the  best  fit 
of  all  locations. 

Aleutian  regional  anomalies 

llie  derived  anomaly  functions  shown  in  Figures  23a  and  25b 
generally  varv  nearly  continuously  with  longitude.  This  result 
nav  appear  to  be  at  variance  with  the  conclusions  drawn  in  previous 
reports  (Chiburis  1968;  Chiburis  and  Aimer,  1970),  whqrein  station 
anomalies,  determined  from  spatially  separated  explosions,  were 
siiown  to  be  essentially  constant,  at  least  across  an  area  approxi¬ 
mately  70  km  x  25  km.  However,  it  must  be  pointed  oat  that  the 
functions  derived  in  this  report  are  made  to  change  continuously 
for  two  reasons:  (1)  most  of  the  events  are  farther  apart  than  the 
distance  over  which  the  anomalies  would  normally  be  expected  to 
be  constant  (probably  about  100  km),  such  that  a  Single  anomaly 
could  not  be  sensibly  used  as  a  correction;  and  (2)  convenience. 

In  spite  of  this,  and  also  to  provide  a  single  value  for  a  finite 
region,  constant  anomalies  within  regions  of  the  Aleutian  Islands 
can  be  determined  from  the  values  shown  in  Figures  25a  and  2ob 
by  averaging  several  anomalies  over  areas  somewhat  larger  than 
desirable.  Because  of  the  lack  of  events,  the  regions  and  the 

o 

-28- 


constant  anomalies  are  arbitrary,  as  the  results  of  determining 
the  step  functions  indicate  in  Figures  25a  and  25b.  Ihe  location 
resuits  obtained  by  using  these  regional  anomalies  are  summarized 
in  Table  XV,  in  which  it  can  be  seen  that  good  clustering  and  low 
standard  deviations  are  being  achieved  to  about  the  sane  degree 
as  in  the  Aleutian  functional  anomaly  case.  Interestingly  enough, 
the  onlv  event  for  which  Network  4  was  not  included  in  tne 
results  is  Andreanof- 1 1 .  This  suggests  that  perhaps  constant 
anomalies  over  fairly  broad  regions  nay  provide  more  stability 
to  networks  with  geometries  like  Network  4. 

In  an  operational  sense,  the  regional  anomalies  thus 
determined  (Figures  25a  and  25b)  are  the  ones  to  apply  when 
locating  events  in  the  Aleutian  Islands  with  any  combination 
of  the  14  stations  in  this  report. 

{  Incidentally,  additional  stations  can  be  added  to  this  14- 
station  network  by  first  locating  one  or  (.better)  several  events 
with  those  stations  having  anomalies.  With  the  resultant  epi¬ 
centers,  an  (average)  anomaly  can  then  be  computed  for  the 
station  to  be  added.  For  subsequent  events  in  that  region,  the 
original  network  and  the  new  station  (or  subsets  of  tnenj  can 
now  be  used  for  locating.  In  this  "bootstrapping"  way,  any 
number  of  new  stations  can  be  incorporated  into  an  existing 
network  such  that  eventually  all  stations  of  any  importance 
would  be  able  to  contribute  equally  and  consistently  to  any 
solution. 

Summarizing  the  results  of  fitting  functions  to  the  obser\ed 
— anomalies  in  various  ways,  we  conclude  that  the  set  of  constant 
regional  functions  derived  are  plausible  and  estimate  the  station 
anomalies  for  events  in  the  Aleutians.  The  anomalies,  when 
applied  as  corrections,  are  believed  capable  of  providing  more 
accurate  relative  locations  for  events  between  the  Fox  Islands 
and  the  .*ear  Islands. 
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SYNTHETIC  EARTH  STUDIES 


All  the  preceding  results  were  obtained  using  real  events 
and  real  arrival  tine  data  receded  by  fourteen  stations.  These 
data  necessarily  contain  tine  errors  which  are  due  to  il  reading 
precision  because  of  less  than  perfect  signal  quality;  2j  unaccount¬ 
able  travel-path  effects  (anomaly  instability!;  and  (5  station 
timing  svstea  errors.  It  would  be  desirable  to  completely  eliminate 
the  last  two  types  of  errors  and  to  be  able,  to  concentrate  on 
investigating  particular  aspects  of  the  functional  anomaly  technique 
in  order  to  verify  or  refute  the  general  method,  lo  acnieve  tins, 
a  data  set  was  synthesized  by  assuming  that  the  Herrin  196S> 
table/  are  perfect ;  from  this  ideal  earth,  theoretical  arrival 
tines  were  computed  at  the  sane  fourteen  stations  for  twelve 
hvpqthetical  events  equispaced  across  the  Aleutians  from  the  Fox 
to  the  Near  Islands,  but  including  the  known  LONG  SHOT  location  in 
the  Rat  Islands.  The  latitudes  and  longitudes  of  the  hypothetical 
events  plus  LONG  SHOT  are  given  in  Table  XVI.  For  simplicity, 
all  events  were  assuned  to  have  surface  foci.  The  effect^  of 
ever-present  reading  error  was  incorporated  in  the  following 
manner:  a  s°et  of  numbers  was  drawn  randomly  from  a  normal  popula¬ 
tion  (Huntsberger,  1961)  by  use  of  random  number  table.  The  set 
was  then  scaled  to  maintain  a  standard  deviation  of  approximately 
0.25  sec  for  each  event,  using  all  fourteen  stations,  by  the 

re  1  a  t  i  on 


Yi  =  '  xj)[0.23(S-l)(sj)  “] 


(attributable  to  readingj 


where  Y?  is  the  random  time  error 
applied’ to  the  ith  station  for  thejth  event,  X;  is  the  nunber 
drawn  fron  the  nornal  population,  X-'  is  the  cean  of  the  X-'s 
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hvent  Parameters 


TABLE  XVI 

for  u  Hypothetical  tpicehtets  Plos  LOSS  SHOT 


Latitude 


Event 

^Decrees) 

Fox  1 

Fox  - 

53. SO  X 

53.00  S 

Fox  3 

s:.4o  X 

\nJ  ■* 

52.00  X 

And  5 

51.60  X 

And  6 

SI. SO  X 

And  ? 

51. SO  K 

Rat  S  (Long  Shot) 

SI. 44  X 

Rat  9 

51.40  X 

Rat  10 

51.60  X 

Near  11 

52.20  X 

Sear  1 - 

S2.80  X 

Longitude 
(Decrees ) 


166. 

00 

b 

168. 

00 

b 

l^O. 

00 

b 

172. 

00 

b 

174. 

00 

b 

176. 

.00 

h 

178. 

.00 

It 

179 

.18 

E 

178 

.00 

E 

176 

.00 

E 

174 

.00 

E 

172 

.00 

E 

\ 


/ 

/ 

/ 

N 

\ 


\ 

\ 


j  2 

for  the  jth  event,  N  is  the  number  of  stations  (14),  and  (S^ 
is  the  variance  of  the  numbers  drawn  for  the  ith  event.  The 
entire  set  of  random  errors  is  given  in  Table  XVII  by  station 
and  by  event.  So  random  errors  were  applied  to  LONG  SHOT, 
because  when  the  anomalies  (determined  with  real  arrival  tines 
are  considered,  random  errors  are  already  included. 

The  location  errors  produced  solely  by  these  reading  errors 
are  shown  in  Figure  26  for  the  sane  network  subsets  as  deiined 
in  the  previous  sections  of  this  report./(Table  Vj.  The  average 
location  errors  for  each  network  are  as  follows: 


Network 

1 
■ y 

5 

4 


Average  Error,  kn 
4.1 


3  .  / 

5.4 

7.2 


These  errors  generally  reflect  a  lower  limit  of  accuracy 
expected  for  each  network  in  locating  events  in  the  Aleutian 
Islands  when  reading  errors  are  on  the  order  of  0.25  sec. 

Svnthetic  functional  anomalies  1 
«  ■■  ■  '  "  " 

Using  the  observed  LONG  SHOT  anomalies  as  constraints, 
piecewise-linear  synthetic  functional  anomalies  were  contrived 
for  each  station  and  are  plotted  in  Figures  2/ a  and  2-b  as 
functions  of  longitude.  The  only  restrictions  imposed  on  the 
synthetic  functions  were  that  they  remain  within  reasonable 
anomalv  ranges  and  that  they  contain  no  large  discontinuities 
when  viewed  as  a  function  of  longitude. 

If  these  anomalies,  referred  to  as  syn-1,  are  assumed  to  be 
caused  by  the  synthetic  earth,  locations  can  be  obtained  using 
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the  four  networks  with  no  corrections  applied.  Symbolically, 
the  arrival  tine  T  at  station  i  for  each  event  is  made  up  of 

the  following: 

\, 


T. 

1 


H  > 
1 


(A.)  , 

v  i'svn-1 


♦  E. 
1 


where  H  is  the  theoretical  tine  from  the  hypothetical  epicenter, 

A  i  is  the  contrived  anomaly  applicable  to  that  event,  and  E 
is  the  random  error  previously  discussed  (c(E-)  =  0.25  sec  for 

'N  '  1 

tfie  set  of  fourteen  station^).  The  location  results  are  shown  in 
Figure  28.  The  average  and  range  of  location  errors  for  the  four 

networks  and  for  fhe  twelve  events  are  summarized  in  Table  XVIII. 

/ 

These  errors  are  absolute  because  the  true  locations  are 
known.  The  wide  range  of  errors  for  a  particular  network  clearly 
shows  the  difficulty  one  can  have  in  attempting  to  understand 
the  anomaly  effect  across  a  region  as  large  as  the  Aleutian 

l _  / : 

Islands,  especially  when  real  events  are  mislocatcd  and  are  at 
different  depths.  Network  3  appears  to  exhibit  the  most  instability 
although  it  has  a  geometry  similar  to  that  of  Network  2.  The 
reason  £s  that  the  stations  in  Network  3  have  the  most  erratic 
and  largest  functional  anomalies. 

\ 

The  contrived  anomalies  and  the  anomaly  variations  for  the 
stations  used  are  not  considered  too  extreme  or  unrealistic. 

Similaf  and  even  larger  variations  have  been  observed  in  the 
station  anomalies  computed  from  the  following  explosion  regions : 
the  Nevada  Test  Site  in  southern  and  central  Nevada,  the  SHOAL 
event  in  western  Nevada,  the  RULISOX  event  in  northwestern  Colorado 
the  GASBUGGY  event  in  northwestern  New  Mexico,  and  the  GNOME  event 
in  southeastern  New  Mexico  (Chiburis  and  Ahner,  1970J.  The  range 
of  anomalies  computed  from  these  explosions  at  some  stations 
exceeds  5  sec,  with  2-3  sec  being  common.  Significantly,  the 


o  o 


■/ 


/ 

v  rf 


distance  between  SHOAL  and  GASBUGGV  is  about  loOO  kr>,  which  is 


of  the  sane  order  as  the  distance  between  the  Fox  islands  and 


Near  Islands 


/ 


The  results  as  shown  in  Figure  28  and  Table  XVIII  were 
obtained  with  four  subsets  of  a  fourteen  station  network.  If 
nore  stations  were  available.  and  other  no  re  stable  network 
subsets  taken,  the  location  of  any  particular  event  in  the 
Aleutians  would  still  be  meaningless ,  except  in  a  general  wav 
(see  the  earlier  section  of  this  report  on  LONG  SHOT  using 
subsets  of  a  529-station  network).  The  point  is  that  if  no 
information  is  available  on  travel-time  anomalies  from  a 
region  of  interest,  tiien  depending  on  the  particular  network 
selected  tor  locating  (even  assuming  the  network. is  reasonablv* 
stable),  the  true  location  errors  ,nay  be  anywhere  from  zero  to 
50  km  or  more.  It  can  be  further  pointed  but  that  in  sole 
regions  the  anomalies  may.  be  distributed  such  that  the  resultant 
errors  would  be  larger.  ^ 

Referring  again  to  Figure  28,  the  results  using  tiie  four 
networks  for  each  event  are  quite  variable,  but  for  any  one 
network  the  shift  vectors  from  event  *to  event  are  co*  durably 
more  consistent;  although  the,  errors  are  large,  the>  vary  slowly 
and  regularly  across  the  region.  For  example,  using  Network  1, 

Fox-1  shifts  about  21  kn  north-northwest  from  thc'Yrue  location, 
and  LONG  SHOT  shi fts^ about  lo  kn,  but  also  to  the  north-northwest. 
Therefore,  Fox-1  ^latlVe^to  LONG  SHOT  is  t>  km  in  error,  etc. 

Table  XIX  gives  the  absolute  and  relative  (to  LONG  SHOT)  errors 
for  all  events  located  by  Network  1.  On  the  average,  the  relative 
errors  afe  reduced  by  a  factor  of  two  over  the  absolute  errors. 

For  Netwprk  3,  which  happens  to  have  the  poorest  distributic-  of 
anomalies  (as  evidenced  by  the  high  average  absolute  error),  the 
same  comparison  relative  (to  LONG  SHOT)  yields  the  results  given  in 
Table  XX. 
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TABLE  XIX 

Absolute  and  Relative  (t<d  LONG  SHOT)  Location  Errors  for 
Hypothetical  Events  Using  Network  1  and  Synthetic  Earth-Model 


Absolute  Location  Relative  Location 


iiv  ent 

Error,  km 

Error,  kra 

I  ox  - 1 

20.8 

S.3 

i  ox-: 

21.9 

6.2 

Fox  -  5 

21.1 

5.6 

And  -  4 

18.3 

2.6 

And  -  S 

14.1 

4.0 

And  -  o 

13.8 

12.9 

And  -  “ 

12.3 

11.7 

LONG  SHOT 

IS. 7 

Rat -9 

22.6 

8.3 

Rat  -  10 

50.7 

15.4 

Near- 1 1 

33.6 

18.6 

Near-i: 

39.9 

24.9 

Ave  rage 

22.1 

10.5 

N 

f 

*/ 


''V. 


. 


y 
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TABLE  XX 


Absolute  and  Relative  (to  LOS'G  SHOT)  Location  Errors  for 
Hypothetical  Events  Using  Network  3  and  Synthetic  Earth-Model 


r* 


Absolute  Location  Relative  Location 


Event 

Error,  kn 

Error,  kn 

Fox-1 

48.5 

21.5 

Fox-2 

42.1 

15.6 

Fox-3 

45.6 

12.2 

And- 4 

40.2 

7 

And- 5 

34.3 

6.1 

And- 6 

26.2 

16.6 

And- 7 

29.1 

15. 5 

LONG  SHOT 

56.9 

'  /)■ 

X  - 

Rat -9 

44.3 

8.2 

Rat-10 

62.8 

25.8 

Near- 11 

67.4 

50.5 

Near- 12 

95.1 

58.1 

Average 

47.7 

19.4 

/ 
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The  re  fore  although  each  network  yields  very  different 
locations  for  an  event  (due  to  different  distributions  of  the 
station  anomalies  within  that  network),  the  average  location 
errors  of  the  set  of  events  relative  to  a  calibration  event  are 
reduced  by  a  factor  of  at  least  two  when  a  constant  network  is 
used.  This  result  further  substantiates  the  conclusions  drawn 
in  a  previous  report  (Chiburis  and  Ahner,  1970)  in  which  the 
relative  location  accuracy  of  2-3  km  for  a  fairly  small  region 
was  shown  to  be  unaffected  by  travel-time  anomalies  if  and  only 
if  a  constant  network  was  used  for  locating  a  set  of  events. 
However,  for  the  synthetic  data  set  used  here,  and  for  the  real 
Aleutian  Islands  data  set  in  general,  it  is  known  that  the 
station  anomalies  across  such  a  large  region  are  not  constant; 
thus  the  relative  accuracy  for  any  network  is  expected  to  be 
lower  than  that  possible  for  the  location  of  events  within  a 
small  region  across  which  the  anomalies  are  constant. 


hquallv  important  as  the  criterion  of  constant  anomalies, 
however,  is  that  of  "network  effect”.  Any  seismic  network  is 
uniquely  defined  by  the  stations  in  it  and  their  geometrical 
relationships  to  the  particular  epicentral  region  of  interest. 
These  relationships  are  functions  of  epicentral  distance  1.  and 
azimuth  0,  and  of  the  travel-time  table,  11(2.),  neglecting  depth 
of  focus.  Specifically,  for  the  i th  station. 


31^(2)  ^s in  9^ 

^  ^cos  ej 

are  the  slopes  of  the  travel-tine/distance  curve  modified  by 
the  sine  or  cosine  of  the  azimuth  depending  respectively  on 
whether  the  rates  of  change  with  longitude  or  latitude  are 

\ 

» 

\ 
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desired.  If  the  region  under  investigation  is  large  enough,  these 
slopes  of  the  travel-time  curve  change  significantly  at  a  station 
sdch  that  the  advantages  of  constant  network  and  constant  anoraalip 
are  considerably  undermined.  This  fact  may  be  demonstrated  with 
the  synthetic  data  by  allowing,  the  known  LONG  SHOT  anomalies  to  ' 
be  the  only  time  errors  across  the  entire  Aleutian  Islands  region. 
In  this  case,  the  arrival  time  data  are  constructed  as 

—  —  } 

Ti  a»i  *  <Ai>LS  •  -  _  j 

{ 

where  (A)LS>i,s  the  LONG  SHOT  anomaly  and  the  other  terms  being 
defined  as  before.  This  means  that  if  the  LONG  SHOT  anomalies  arte 
applied  as  corrections  to  the  arrival  time  data,  all  location  / 
errors  obtained  with  any  network,  constant  or  otherwise,  would  i 
be  zero.  If,  on  the  other  hand,  constant  networks  are  used  but' 
no  anomalies  are  applied,  the  results  are  as  shown  in  Figure  29. 
Here  it  can  be  seen  that  Networks  1  and  2  have  only  slight  network 

CC  -  ’’Yl'*’ 

effects  on-relative  location  accuracy  due  to  slope  perturbations. 
Network  3,  however,  has  a  strong  network  effect  in  that  the  errors 
have  similar  directions,  but  the  magnitudes  differ  by  a  factor  of 
almost  3  between  the  Fox  and  Near  Islands.  The  effect  on  Network  4 
produces  errors  which  change  direction  by  about  60°  and  magnitude 
by  a  factor  of  2/3  between  the  Fox  and  Near  Island.  These  shift 
differences  are  due  to  network  geometry,  to  the  network  relation¬ 
ship  to  the  Aleutian  Islands  region,  and  to  the  non-linearity  of 
the  travel- time  curves,  because  the  network  itself  and  the  station 
anomalies  are  held  constant;  the  only  parameters  which  are  variable 
are  the  station  distances  and  azimuths  to  the  twelve  events.  If 
constant  anomalies  were  the  only  criterion  necessary  to  maintain  \ 
relative  accuracy  over  a  broad  region,  the  location  shifts  for  a 
constant  network  would  all  have  had  the  same  magnitude  and 
direction.  This  is  true  over  a  small  portion  of  the  Aleutian 


f 


i 

region,  in  which  the  location  shifts  an ’  . 

it  is  not  true  in  general  “to  cvcnt  g,  which  is 

using  Network  3,  And-' '  a"  3  kn  rcspect i ve ly .  whereas 

LOSC  SHOT.  H-;-  s;fr^ire  to  event  S  have  errors  of  10 
events  Fox-1  and  we.  -  arc  far  larger  than 

and  16  kn  respectively.  The  lat  jmU  csough  that 

the  5-5  kn  errors  expected  the  g  .  changc. 

the  travel- tine /distance  slopes  do  or  si  ni  ^  ^ 

Therefore  if  *  ““‘“‘..“‘“gio^but  are  not  applied, 
are  constant  throughout  B  *  f  traVcl-ti=e  slopes 

relative  accuracy  dininishes  as  a  to  the 

•  „  distance  between  events  due 

with  increasing  <  *  anonalies  are  applied  as 

network  effect.  Of  course  t  those  due 

corrections ,  -t-ouU^no  location^  ^  ^  ^  ^ 

to  reading  Pre‘151°  *  data  which  inciude  the  synthetic 
Table  XIX  and  XX,  usi  g  kt  least  sone  inprovenent 

functional  anonalies.  one-  ^  the 

in  relative  accuracv.  .reater  than  that  produced 

anonalies  are  producing  an  effect  .re 

by  network  geonetry.  anonalies  are  included 

When  the  synthetic  function  •  assuned  that  the  LONG  SHOT 

.  if  it  is  erroneous  1>  assume 

l^iin^-id  across  the  entire  -- fans  Is  lands  region . 

—  -  r^-elrks.  The  arrival  tine 

locations  can  be  nade 

data  now  are 


Ti  -  ^Ai^svn-1 


♦  F. 


(-Vls 


.  with  this  relation  can  also 

the  process  of  “^^/“t^lv^tracting  vectorially  the 
e  graphically  obtained  b>  smpl.  * 
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shifts  shown  in  Figure  29  (LONG  SHOT  anonalies)  fron  the  shifts 
in  Figure  2$  i.syn-1  anoaalies)  for  each  event  and  for  each 
network.  For  Fox- 1  and  for  Network  4,  this  yector  process  is 
shown  in  Figure  50).  The  average  location  results  for  each  net- 
worif  are  given  in  Table  XXI,  which  also  includes  for  conpanson 
the  averages  already  given  in  Table  XVIII  obtained  when  no 
anomalies  were  applied. 

For  all  networks,  the  average  errors  are  reduced,  and  except 
for  Network  2,  which  has  a  low  average  to  begin  with,  the  reduc- 
tiohs  are  significant.  These  results  suggest  that,  although  the 
LONG  SHOT  anonalies  are  not  completely  valid  across  the  entire 
region,  they  are  at  least  partially  effective  in  renoving  sone 
of  the  gross  bias  errors.  Stated  differently,  the  effect  of  the 
LONG  SHOT  anonalies  is  greater  than  the  effect  of  the  synthetic 
anomaly  variations;  when  the  LONG  SHOT  anoaalies  are  applied  as 
corrections,  first-order  bias  effects  are  eliminated  but  second- 
order  bias  effects  renain.  These  second-order  errors  renain  even 
for  those  events  adjacent  to  LONG  SHOT  (And- 6,  And-/,  Rat -9, 
and  Rat-10),  where  the  anonalies  would  be  expected  to  be  approxi- 
nately  valid,  since  the  synthetic  anoaalies  were  Bade  to  change 
only  slowly  as  a  function  of  distance  away  fron  LONG  SHOT. 

Figure  31  shows  the  results  of  applying  the  LONG  SHOT  anomalies 
to  the  synthetic  data  and  locating  with  the  four  networks. 

Table  XXII  sunnarines  the  complete  results.  The  errors  in  the 
mean  location  (snail  circles  on  the  Figure)  are  less  than  10  kn 
for  events  4  through  9;  for  events  6  through  11  the  clustering 
errors  are  reduced  to  less  than  6  kn  and  the  standard  deviations 
are  acceptably  low.  The  errors  in  the  nean  location  indicate 
absolute  accuracy,  while  clustering  errors  indicate  relative 
accuracy.  Therefore,  as  noted  earlier,  clustering  and  low- 
standard  deviation,  alone  or  together,  are  necessary  but  not 
significant  to  guarantee  improved  absolute  accuracy. 
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TABLE  XXI 

Average  Location  Errors  Using  LONG  SHOT  Anomalies 
with  Synthetic  Earth-Model 


\ 


Network 

\ 

Average  Error  km. 
Using  LS  Anomalies 

Error  Range,  kn 

Previous  Average 
Error,  km,  Using 
no  Anomalies 

1 

11.4 

'  3.7  - 

-  23.8 

22.1 

2 

10.6 

1.4  - 

■  16.9 

12.5 

3 

20.0 

9.7  - 

■  42.8 

47.7 

4 

21.4 

5.2  - 

■  49.1 

32.6 
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TABLE  XXII 

Average  Location  Shifts  of  Mean  Epicenter,  Standard 
Deviations  (c),  and  Clustering , Ef feet  Using  LONG  SHOT 
Anomalies  with  Synthetic  Earth-Model 


is 


Error  in 

Average  Shift 

Mean 

Event 

Mean  Location 

Fron  Mean 

a 

Fox  1 

12.7 

25.3 

0.85 

* 

Fox  : 

13.0 

13.9 

0.87 

Fox  3 

13.1 

18.0 

0.81 

\nd  4 

6.8 

10.7 

0.64 

And  5 

1.0 

10.2 

0.54 

And  6 

10.0 
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Figure  30.  Vector  addition  for  resultant  location  shift. 
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Furthermore,  the  LONG  SHOT  anomalies  appear,  to  be  valid 
only  for  the  events^  immediately  adjacent  to  it  (events  7  and  9 
at  distances  of  about  180  km  and  80  km  respectively  from  LONG 
SHOT).  The  mean  location  error  of  event  7  is  about  10  kn;  that 
of  event  9  about  5  kn.  This  result  is  in  general  agreement  with 
previous  conclusions  concerning  the  size  of  an  area  across  which 
a  constant  anomaly  is  good  enough  for  accurate  location  (Chiburis 
and  Ahner,  1970). 

Anomaly  interpolation  using  synthetic  data  N 

The  location  of  LONG  SHOT  is  known.  If  some  location  is 
assumed  for  event  1  (Tox  Islands)  and  for  event  12  (Near  Islandsj  , 
and  if  the  anomalies  are  computed  from  the  three  locations  with 
linear  fuhctions  drawn  between  them,  can  improvement  be  gained  in 
accuracy  and  in  standard  deviation  over  that  obtainable  if  no 
anomalies  are  applied  to  the  intervening  events?  Restated,  can 
one  combine  the  anomalies  from  a  true  location  and  from  wrong 
locations  and  still  expect  to  improve  location  capability?  In 
the  following  discussion  we  will  answer  these  questions  in  the 
affirmative. 

Using  the  synthetic  earth  in  which  the  anomalies  are  known, 
Fyx-1  is  deliberately  dislocated  by  0.2°  north  and  Near-IE  by 
0.1°  south.  Computing  the  anomalies  fr'bra  these  two  fiducial 
locations  and  from  the  known  LONG  SHOT  location,  bilinear 
functions  are  fitted  to  the  three  constraints,  as  are  shown 
in  Figures  32a  and  32b.  Compare  these  functions  with  the  synthetic 
functions  originally  contrived  (Figure  27a  and  27b). 

Locating  the  twelve  hypothetical  events  with  the  four  net¬ 
work  ^subsets  and  the  bilinear  functions  (randoci  reading  errors 
are  also  included),  gives  the  results  shown  in  Figure  33.  Because 
events  1  and  12  were  deliberately  given  a  location  bias  of  0.2°N 
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3nd  0.1°S  respectively,  in  order  to  discuss  accuracies  across 
the  region,  the  other  events  need  to  be  adjusted  by  an  amount 
proportional  to  their  locations  relative  to  the  calibration 
events.  The  adjustment  relation  is  shown  in  Figure  34.  The 
adjusted  fiducial  epicenters  are  shown  as  circles  in  the  pre¬ 
vious  Figure  33. 

■) 

The  location  errors  relative  to  the  adjusted  epicenters 
and  standard  deviations  are  given  in  Table  XXIII.  Compared  to 
the  no-anonalv  results  with  the  same  synthetic  data  (Table  XVIII), 
significant  improvement  is  obtained  in  relative  accuracy  by 
assuming  the  locations  of  "unknown"  events  and  tying  the  anomalies 
to  those  of  known  explosions.  The  reason  that  the  events  And-6 
and  And-7  do  not  locate  well  is  due  to  the  nature  of  synthetic 
functions  derived;  many  of  the  piecewise- 1 inear  functions  were 
made  to  change  sharply  between  events  5  and  8,  which  results  in 
poor  estimates  of  the  anomaly  when  assuming  strict  linearity 
between  the  calibration  events.  Therefore,  across  a  region  the 
sice  of  the  Aleutian  Islands,  the  technique  of  anomaly  inter¬ 
polation  is  valid  and  will  yield  improved  relative  locations. 

This  conclusion  bears  directly  on  the  interpolation  results 
obtained  with  real  data;  viz,  when  several  events  are  assumed 
to  be  located  correctly  and  linear  anomalies  are  fitted  between 
them,  the  results  are  consistent  both  from  clustering  and 
standard  deviation  considerations  (Figure  24  and  Table  XIII). 

Step-off  technique 

A  brief  discussion  is  now  given  of  a  method,  referred  to  as 
the  "step-off"  technique,  that  attempts  to  increase  the  area 
over  which  the  known  LONG  SHOT  anomalies  remain  valid.  The  obser¬ 
vations  which  suggested  the  possible  applicability  of  this  tech¬ 
nique  were  that  the  distribution  of  -'zero-mean  tine  errors  resulting 
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fron  the  least-squares  solution  followed  fairly  closely  the 
distribution  of  functional  anonalies  contrived  for  the  Aleutian 
Islands  region.  For  exanple,  using  Network  1  and  the  synthetic 
data,  the  set  of  12  hypothetical  epicenters  was  located  using 
onlv  the  LONG  SHOT  anonalies  as  corrections  across  the  whole 
region;  the  residuals  for  stations  HN-ME,  STUf  and  TrO  resulting 
fron  the  least  squares  solutions,  and  the  anonalies  originally 
contrived  for  these  stations  are  plotted  as  functions  of  longitude 
m  Figure  35.  As  can  be  seen  in  the  Figure,  these  residuals 
estinate  quite  well  the  anonalies  actually  incorporated  into  the 
arrival  tine  data,  especially  in  the  vicinity  of  LONG  SHOT.  At 
larger  distances,  the  residuals  begin  to  diverge  fron  the  contrived 
anonalies  as  expected.  Therefore,  if  one  applies  the  LONG  SHOT 
anonalies  to  an  adjacent  event  (say  event  Aj  and  then  attributes 
the  cause  of  the  residuals  resulting  fron  the  location  of  event  A 
to  a  real  change  in  the  anonaly  between  the  two  events,  these 
errors  can  be  algebraically  added  to  the  LONG  SHOT  anonalies, 
and  the  sun  applied  as  an  anonaly  correction  to  the  next  event 
(say  event  B) .  In  this  way,  one  can  "step-off"  away  fron  LONG  SHOT, 
event  by  event,  and  hopefully  inprove  the  bias  nornally  observed 
between  the  Fox  and  Near  Islands.  The  i’th  station's  anonaly  to 
be  applied  to  the  first  event  in  distance  aWay  fron  LONG  SHOT  is 
just  the  LONG  SHOT  anonaly: 
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These  anonalies  yield  a  location  of  event  A  with  a  set  of  least- 
squares  errors  .  The  step-off  anonaly  to  be  applied  to  the 

second  event  away  fron  LONG  SHOT  is 
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Although  the  application  of  the  step-off  anonalies  results 
m  a  considerably  reduced  standard  deviation  in  the  solutions, 
all  of  the  final  locations  were  found  to  be  nearly  identical 
within  1/2 'kn)  to  those  obtained  by  applying  only  the  LONG  SHOT 
anonalies.  In  other  words,  applying  a  set  of  residuals, 
which  results  fron  the  location. of  one  event,  to  the  set  of 

arrival  tines  of  an  adiacent  event  has  a  negligible  effect  on 

“  *  *  ' 

the  location,  ever  though  these  applied  errors  are  large.  Table 
XX i \  gives  the  least-squares  errors  from  the  final  locations 
obtained  by  using  the  LONG  SHOT  anomalies;  the  three  standard 
deviations,  =Q,  rLS,  cLS+£,  are  respectively  those  obtained 
with  no  anona'ly  corrections,  those  with  the  LONG  SHOT  anomalies, 
and  those  with  the  LONG  SHOT  anomalies  plus  the  least-squares 

j  • 

errors  fron  the  adjacent  event.  The  least-squares  errors  in 
locating  event  7  (Andreanof)  with  the  LONG  SHOT  (event  8) 
anonalies  are  given  in  the  first  column  of  Table  XXIV.  This 
set  has  a  standard  deviation  of  0.201  sec.  (The  standard 
deviation  obtained  by  applying  no  anomalies  is  shown  as  eQ.) 
Adding  these  errors  to  the  LONG  SHOT  anomalies  and  locating 
event  6  results  in  a  standard  deviation  of  0.408  sec  instead 
of  0.4'S  sec  as  obtained  with  the  LONG  SHOT  anonalies  alone. 

When  locating  event  1  (Fox)  using  only  the  LONG  SHOT  anonalies, 
the  standard  deviation  is  1.286  sec;  if  the  least-squares  errors 
fron  event  2  are  stepped-off,  the  standard  deviation  is  reduced 
to  0.414  sec.  For  the  synthetic  data  set,  the  standard  devia¬ 
tion  when  the  LONG  SHOT  anonalies  are  applied,  :?LS,  compared 
to  the  standard  deviation  without  anomalies,  c Q,  is 


TABLE  XXIV 

Set  of  Least-Squares  Errors  Applied  in  Step-Off 
Technique  and  Resultant  Standard  Deviations.  (See  Textj 

\  V 

LEAST-SQUARES  ERRORS  (IN  SECONDSj 


Station 

AXD7 

AXD6 

AND  5 

COL 

-.23 

-.54 

-.84 

XP-XT 

-.10 

.22 

.64 

XOR 

.47 

.90 

1.12 

LA010 

.11 

.18 

.08 

UBOl 0  , 

-.27 

QO 

fN 

• 

1 

-.61 

RK-OX 

-.13 

.54 

.53 

TF060 

.02 

-.32 

.25 

TUC 

.19 

.76 

1.11 

NM006 

.17 

-.50 

-.71 

XUR 

-.02 

-.46 

.01 

HX-ME 

-.01 

-.15 

-.77 

XDI 

-.01 

.26 

.08 

QUE 

-.25 

-.30 

-.60 

STU 

.08 

-.31 

-.28 

oo 

.941 

1.074 

1.2^4 

oLS 

.201 

.478 

.672 

oLS*E 

.201 

.408 

.357 

AND  4 

FOX  3 

FOX  2 

FOX  1 

-1.01 

-1.58 

-1.51 

-1.81 

.63 

•  .75 

.64 

1.25 

1.61 

2.31 

2.46 

2.23 

.89 

.81 

1.08 

.98 

-.23 

-.38 

.24 

-.56 

.25 

.74 

.55 

.72 

.17 

.14 

-.39 

.30 

.57 

.45 

.45 

.28 

-.82 

-.78 

-.99 

-.93 

-.13 

-.31 

-.21 

.23 

-.84 

-.70 

-.58 

-.37 

.26 

.75 

1.09 

1.16 

-.69 

-1.74 

-1.96 

-2.25 

-.65 

-.46 

-.85 

-1.34 

1 . 366 

1.470 

1.536 

1.666 

.766 

1.072 

1.168 

1.286 

.357 

.444 

.304 

.414 

reduced  signif icantlv-in  the  vicinity  of  the  LONG  SHOT  explosion, 
which  indicates  that  the  LONG  SilOT  anomalies  are  estimating  well 
the  contrived  anomalies;  the  improved  locations  for  these  events 
l  iable  XXI  I )  bear  this  cut.  But  all  of  the  standard  deviations 
with  the  step-off  errors  applied  compared  to  the  are 

also  reduced  but  the  location  errors  are  not.  This  would  indicate 
onlv  tli at  the  contrived  anomaly  changes  are  being  estimated  well 
bv  the  step-off  process,  but  since  the  location  errors  remain 
identical  to  those  obtained  with  the  LONG  SHOT  anomalies,  nothing 
is  being  gained  except  artificially  low  standard  deviations. 

The  preceding  exercise  was  made  merely  to  demonstrate  the 
futility  of  using  a  set  of  computed  errors  (residuals)  from  a 
network  solution  in  an  attempt  to  improve  the  location  of  a 
subsequent  and  nearby  event.  The  reason  this  cannot  work  is  that 
since  the  set  of  errors  has  a.  minimum  value  in  a  least-squares 
sense,  and  since  the  network  is  uniquely  defined  by  the  stations 
within  it,  the  set  of  errors  applied  to  some  event  already 
located  irv  the  region  by  that  network  has  a  mathematically  null 
effect  as  fat  as  further  perturbations  are  concerned.  Stepping-off 
slowly  in  this  way  across  a  broad  region  keeps  this  effect 

negligible,  and  the  resultant  locations  remain  the  same  as  if 

*  . 

no  corrections  wpre  applied. 
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CONCLUSIONS 


This  study  deals  with  travel-tine  anomalies,,  networks, 
and  location  consistencies  fron  two  viewpoints  real  data 
and  hypothetical  data  for  Aleutian  Islands  everts.  The  hypo- 
thethical  data  were  used  to  validate  sone  of  the  techniques 
enployed  on  the  real  data.  Physically,  the  technique  of 
functionalizing  and  regionalizing  the  anomalies  works  well 
enough  to  justify  further  study  and  general  applicability. 

No  implication  is  made  as  to  the  ease  with  which  it  can  be 
accomplished,  but  forefaost  anong  the  criteria  for  using  it 

i  • 

is  that  the  raw  data  set  nust  be  high  quality;  that  is,  the 
arrival  tines  must  read  as  unambiguously  as  possible.  This 

i 

usually  involves  a  great  deal  of  intrastation  checking  and 
interstation  correlation.  The  earthquake  signals  from  Aleutian 
Islands  events  dre  frequently  subtly  complicated:  double  events 
separated  by  2r5  seconds,  usually  very  different  in  magnitude; 
low  amplitude  first  notions  (which  may  be  called  precursors) 
not  normally  observed  at  noisier  or  low-iaaiini f ication  stations, 
etc.  Since  studies  of  anomalies  and  location  techniques  should 
not  concern  themselves  with  the  ability^of  a  seismic  analyst 
to  read  seisraograms  correctly , the  majoA  effort  of  future  studies 
must  be  to  obtain  as  clean  a  data  set  as  possible  using  all 
available  techniques. 

In  using  the  real  data  for  Aleutian  Islands  events,  the 
complications  introduced  in  the  location  stabilities  are 
undoubtedly  due  in  part  to  a  depth  effect.  If  the  anomalies 
are  caused  bv  the  integration  of  small  imperfections  in  the 
earth  model  along  a  unique  travel  path  (which  is  then  related 
to  a  region/station  description),  there  is  good  reason  to 
expect  the  anomalies  to  be  dissimilar  when  computed  from  two 
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events  m  the  same  region  but  at  significantly  different  depths. 

It  is  not  presently  known  what  sice  of  depth  differences  are 
to  he  considered  significant,  but  the  significance  would  probable 
depend  on  the  particular  stations  used.  However,  before  sensible 
studies  can  be  made  relevant  to  this  depth  effect,  all  possible 
efforts  must  be  made  to  eliminate  (1)  arrival-time  errors; 

(_ j  i  epicenter  errors;  (5)  depth  errors;  and  (4)  network  insta¬ 
bilities.  To  eliminate  the  first  of  these  requires  labor  and 
patience.  To  eliminate  the  second  may  be  impossible,  except  in 
the  case  of  explosions  or  of  using  well-distributed  local 
networks  (<100  km  epicentral  distance).  To  eliminate  the  third 
requires  the  success  of  the  first  two  plus  the  clear  observation 
of  pP.  The  fourth  may  be  el iminated , . or  certainly  be  made  manage¬ 
able,  by  selection  of  suitable  networks. 

Neglecting  d^th  effects,  the  analysis  of  the  real  and 
hypothetical  data  using  various  interpolation  and  adjustment 
schemes  permits  an  estimate  of  resultant  relative  location  error 
to  be  placed  at  5-10  km  throughout  the  Aleutian  Islands  for  the 
seventeen  events  studied.  The  absolute  location  error  is  probably 
about  the  same,  because  LONG  SHOT  provides  bias  control  in  the 
vicinity  of  the  Hat  Islands  and  because  all  functional  and 
regional  anomalies  were  tied  to  LONG  SHOT. 

The  regional  anomalies  finally  determined  for  the  Aleutian 
Islands  yield  consistent  location  patterns  and  least-squares  time 
fits  using  four  subsets  of  a  fourteen  station  network. 

[hat  epicenters  may  be  mislocated  in  the  Aleutian  Islands 
regions  was  demo. is trated  by  taking  various  stable  subsets  of  a 
519-station  network  and  the  LONG  SHOT  explosion.  The  .variable 
locations  obtained,  any  of  which  could  have  been  reported  had 
LONG  SHOT  been  an  earthquake,  produce  a  computed  anomaly  range  at 
the  teleseismic  station  NP-NT  of  more  than  11.4  sec.  This  result 


implicates  epicenter  mislocation  as  the  chief  cause  of  obserrei 
anomaly  scatter. 

Finally  we  have  shown  that  arrival  times  in  published 
bulletins  are  frequently  in  error  by  1  second  and  more.  The 
causes  of  these  errors  are  probably:  (1)  m .ssing  first  i  ti'n 
in  complicated  signals  (such  as  double  events),  and  C)  «».  . 

first  motions  (as  at  noisy  stations).  Any  future  studies  o. 
anomalies,  networks,  location  accuracies,  and  earth  ode  s 
must  begin  with  an  analysis  of  raw  seismograms  in  order  to 
place  any  significance  on  the  redation  of  anomalies  or  anomaly 

variations  to  physical  processes. 
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